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When an informationally incomplete set of observables is considered there are several solutions to 
' the quantum state reconstruction problem using von Neumann measurements. The set of solutions 

M ' are known as Pauli partners, which are not easy to find even numerically. We present, in a self- 

contained paper, a new way to find this solutions using the physical imposition operator. We show 
that every Pauli partner is an attractive fixed point of this operator, which means that we can 
find complete sets of Pauli partners very efficiently. As a particular case, we found numerically 24 
mutually unbiased bases in dimension A = 23 in less than 30 seconds in a standard PC. We hope that 
the algorithm presented can be adapted to construct MU Constellations, SIC-POVMs, Equiangular 
Tight Frames and Quantum t-Designs, which could open new possibilities to find numerical solutions 



INTRODUCTION 



^ ■ to these open problems related with quantum information theory. 

^H. Keywords; Quantum state reconstruction. Mutually unbiased bases, Pauli partners. 

PACS: 03.67.-a;03.67.Ac 

> 
m 

(N 
rn 

rn 

. Quantum state reconstruction from complete sets of physical distributions is an important open 

problem in the foundations of quantum mechanics. It has been proved that it belongs to the 
complex category NP-(Hard)[l|, where NP means non polynomial. That is, this is a problem 
' harder than NP-problems. An algorithm is NP when the time required to solve it increases faster 

I than any polynomial function depending on the relevant parameters of the system, for example, 

the space dimension. Sometimes, there exists polynomial algorithms that can solve NP problems 
efficiently Q. In this case, we say the problem is NP-EASY. For example, the sort operation is a 
NP problem, but there is a polynomial operation (quick sort), that is able to solve it in polynomial 
time. A review of NP problems that can be reduced to polynomial ones can be found in [3]. In our 
work, we could not find a polynomial algorithm for quantum state reconstruction. However, we 
have an exponential algorithm that has a very small exponential factor, what means that we can 
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reconstruct efSciently quantum states, and complete sets of Pauli partners, in high dimensional 
Hilbert spaces. For example, we found numerically 24 mutually unbiased bases in a 23-dimensional 
Hilbert space in 27 seconds with a standard PC. Here, we present a new way for studying quantum 
state reconstruction using dynamical system theory. We show that different states containing 
the same information about given observables, that is Pauli partners, are attractive fixed points 
of the physical imposition operator. We will present an algorithm able to find complete sets of 
Pauli partners when informationally incomplete 37[ sets of observables are considered and a unique 



solution when the observables are informationally complete. As a particular case, our algorithm is 
able to recognize when a set of observables is informationally complete. 

This paper is organized in the following way: in Section II, we define the quantum state re- 
construction problem and we mention some related open problems. In Section III, we present an 
algorithm for state reconstruction and we prove several properties using dynamical systems theory. 
Additionally, we give a sufhcient condition for defining informationally complete set of observables. 
In Section IV, we discuss numerical results and we mention some possible applications for the al- 
gorithm, and in Section V we resume our work. In the appendices there are some proofs an we 
define a generalization of Hellinger metric. 



II. PAULI PROBLEM 

The pure state of a quantum system is a ray in the Hilbert space, and it contains the complete 
information about every relevant physical property of the system. In this work, we consider the 
quantum state reconstruction from von Neumann measurements, that is, from complete sets of 
eigenvalues distributions corresponding to several observables. 

In the laboratory, we are able to make measurements that contain information about the quan- 
tum state. However, these measurements are not the projections of the state in a basis of the 
Hilbert space, like happens in classical mechanics. What we are able to measure are only the 
amplitudes of these projections, and this means that we should solve a set of non linear equations 
in order to find the state of the system. When we consider distributions from an incomplete set 
of observables, these equations do not have a unique solution, and this situation is precisely the 
Pauli problem: there are more than one state with the same partial information about the phys- 
ical system, and the problem is to find the complete set of them. This situation was suggested 
originally by W. Pauli in a footnote]^, and the problem is also known as the quantum state re- 
construction problem. Different states containing the same set of eigenvalues distributions of the 
quantum system are named Pauli partners. The original version of Pauli considers position and 
momentum observables, and there are several examples of different quantum states with the same 
distributions for these observables 0-Q]- We particularly remark Theorem 2.5 in Q| where it is 
proved that there exists wave functions with an infinite number of Pauli partners. Today, we know 
that the Pauli problem appears when we do not consider an informationally complete P] set of 
observables. In the case of position and momentum, we know that their eigenvectors basis are not 
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the maximal set of mutually unbiased bases (MUBs). For example, in a iV-dimensional Hilbert 
space we need — 1 observables beside position and momentum [10] to reconstruct every quantum 
state efficiently ill]. In quantum information theory, position and momentum bases are known 
as the eigenvectors bases of the shift {X) and phase (Z) operators. Shift and phase operators 
are physically relevant, because they are the generators of the translation and impulsion of states 
in a unity of magnitude. Also, these operators are very close to two important open problems in 
quantum information theory: existence of MUBs and SIC-POVMs in arbitrary dimensional Hilbert 
space. 

Let us give an introduction to the quantum state reconstruction problem. Let A, B be two 
observables and {(pk\k=i..N i{4'p\p=1--N be their eigenvectors basis. The information about A and 
B contained in the pure state $ is given by the eigenvalues probability distributions 

pt = \{ipk.n\k^l..N (1) 

and 

-|(0p,<i>)|2,p = l..iV (2) 

respectively. The state reconstruction problem for the observables A and B consists in finding $ 
from the knowledge of the distributions {p^} and {p^}- 

We remark that Pauli partners exist in every finite dimensional Hilbert spaces, even in the 
smallest dimension N = 2. For example, a 2-dimensional state $ have the same position and 
momentum distributions than $* ^14]. The spin s = 1/2 version of this problem do not have Pauli 
partners because the observables Sx, Sy, Sz are a maximal set of mutually unbiased operators, that 
is, observables with mutually unbiased eigenvector basis. Pauli partners disappear because these 
eigenvectors basis determine a complete set of Mutually Unbiased Bases (MUBs) in dimension N = 
2, and it has been proved that complete set of MUBs are optimal for quantum state reconstruction 



11| . However, for s 7^ 1/2 the spin eigenvectors basis are not a set of MUBs and there exists a null 



17[, we found analytically 



measure set of states having Pauli partners Q, [l^ ■ In a previous work 
the complete set of states with Pauli partners corresponding to a particle with spin s = 1 predicted 
by Amiet-Weigert in jl^ using our physical imposition operator. 

Another interesting problem related to quantum state reconstruction is to find the maximal 
set of MUBs. Two orthonormal bases defined in a N dimensional Hilbert space {(pk}k=i..N and 
{4>p}p=i..N are MUBs if 

|(^fc,<^p)|2 = l,VA:,p=l..Ar. (3) 

It is well known that there exists a maximal set of iV + 1 MUBs when iV is a power of a prime 
number and at most + 1 MUBs in every dimension [10], but the maximal set of MUBs in non 



power of prime dimensions is still unknown, even in the lowest dimensional case TV = 6 181-1201] . The 
MUBs problem is a particular case of quantum state reconstruction problem, that appears when 
we consider MUBs eigenvectors basis for the observables and when every probability distribution 
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considered is flat. A distribution is named fiat when every probability take the same value 
For this reason, our algorithm may be applied to try solve the MUBs problem. 



III. ALGORITHM FOR STATE RECONSTRUCTION 



In finite dimensions, the quantum state reconstruction problem is defined as follows: Suppose 
that we have a quantum system with m incompatible observables. Each one has an eigenvectors 
basis {"ySfclfc^i "™, where N is the dimension of the Hilbcrt space. Suppose also that we have 
realized measurements of the m observables on an ensamble prepared in the same unknown state 
$, obtaining the eigenvalues distributions {pi.}, r = I..m, k = 1..N. Then, the Pauli problem 
implies to find the complete set of solutions for $ satisfying the set of non linear equations 

pl^Mk,^)\',yk = l..N,r^l..m. (4) 

The maximal set of different states satisfying Eq. (jH) determines the complete set of Pauli partners 
of the system. Its general solution is still unknown and only a few particular cases were solved or 



partially solved, for example 13,Q, [l^Hl- Next, we will define a non linear operator which will 



be an useful tool for finding complete sets of Pauli partners. 



A. Physical imposition operator 



The physical imposition operator is a natural way for transmitting to any state the complete 
eigenvalues distribution obtained from measurements over a system. To understand how it works 
let us analyse the increase of knowledge that has an observer about an unknown ensemble when 
he realize measurements. We should do not confuse this situation with the knowledge that has 
the observer about the collapsed state, but with the partial knowledge acquired about the original 
state, unperturbed by measurements. A way to quantify the partial knowledge about an ensemble 
is to consider the volume of all possible states that are compatible with the information collected 
about the state. Before any measure, the observer has no knowledge about the system. So, he is 
not able to rule out any Hilbert space element. In this situation, his best option is to choose a 
state at random, namely ^Pq. Obviously, this state has any property about the system. In order 
to obtain a better description of the quantum system we must impose some information on 
that should be obtained from measurements. Then, let us suppose that the observer considers an 
observable A with eigenvectors basis {(pk}k=i..N and he realizes a standard tomography process 
over an ensemble prepared in the unknown state $. After measurements, the information obtained 
is the distribution 

= Vfc = l..iV. (5) 

In the measurement process, we obtain information about the system and, after that, all states in 
the Hilbert space are not equally probable for describing the original state. Suppose that we want 
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to establish a connection between the state chosen random, that contains nuU information about 
the system and the state ^Pi, having the complete information about the probability distribution 
{pk} and no more information than this about the system. Of course, there are infinite options for 
^1 but, a simple one involves the following operator 

DEFINITION III.l Let A be an observable with eigenvectors basis {(pk}k=i..N cind <!', 5*0 G 
Then, we define the Physical Imposition Operator as 

T^**o^f]|(^,,$)|i^^^,. (6) 

Then, our state is given by Ta^^'q. This operator is well defined for every quantum states 
except when = y^k, for any k — 1..N. When this happen, we replace ||^^'^°|| by the unity. Let 
us analyze in three steps how T/i$ acts on a state 

1- Expands "^o in the eigenvectors basis {(fik} '■ X]fcLi(V'fc, '^o}fk- 

2- Removes the information that ^I^o contains about A: X^fcLi ^fc- 

3- Imposes the information about A obtained in the laboratory: X^fcLi lifki $)! 'Pfc- 

Notice that in the step 2 we have not a quantum state, because the vector is not normalized. 
However, in the step 3 normalization is restored because the distribution imposed is normalized. 
The physical imposition operator is non linear, idempotent and preserve norms. These properties 
can be proved very easily from the Definition IIII.ll However, it is neither a projector nor an 
unitary operator due to the nonlinearity. The physical imposition operator gives us the first 
approach to the state <i> of the system. Notice that the complex unitary phases in the state 'i'o 
remind unchanged when Ta^ is applied. The reconstruction problem would be solved if we could 
define an observable B that let us find the unknown phases in the laboratory but, as we know, it 
is not possible in general. This is the main reason why the quantum state reconstruction problem 
exists. The knowledge about only one observable is not enough for the reconstruction of the state, 
except in the trivial case when the system is prepared in an eigenvector state of the observable. 
Then, we needed to take into account additional information of, at least, a second observable B 
that does not commute with A with the aim to fix the state. In order to do this we consider a 
second physical imposition operator Tb$, related to the observable B, given by 

T..*i = El(^,'J>>l^^^. (7) 

where TTp = \{4>p, is the eigenvalues distribution of B in the state <E> and 4>p are the eigenvectors 
of B. The state ^'i = Ta^^'q has more information than ^I'o about the system and the next step 
is natural, we should consider the state ^1*2 ~ Tb^^E'i = Tb^Ta^^o- Notice that ^['2 contains the 
complete information about the distribution {TTp} but a partial information about the distribution 
{pk}, because the imposition operator Tg^ destroys the information about A when the modulus of 
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the coefficients of 4'i in tlie basis {(/>p} are replaced. However, part of tlie information about {pk} 
remains in the unchanged phases when is applied. It is important to remark that the single 
imposition operators are idempotent, because they exhaust the information about an observable, 
but the multiple imposition operator = Tbc^Ta^ does not exhaust the information about A 

and B, and then Tab* is not idempotent. So, we can define the sequence = (Tas$)"^o and 
hope that the successive impositions could give a state with the complete information about A and 
B. If we have more than two observables we should consider the sequence = (rA_BC..*)"^o- 
In every step of the sequence the unknown phases gain information about the set of distributions 
and, intuitively, we hope that after an infinite number of steps this insistent process converges to 
a solution. In order to study the convergence of this sequence we must define a metric. Basically, 
we need two metrics because 1) We need to know when a sequence converges, and 2) We need to 
know when two states are partners. Therefore we need a metric in the space of Hilbert space rays 
and another metric in the space of eigenvalues distributions. In order to analyze the convergence 
of a sequence of quantum states we cannot take into account the usual metric in Hilbert space, 
because quantum states are defined up to a global complex phase. The usual distance between the 
states $ and e*"vl' is given by (5($, e'"^') and it depends on a. A good choice is the Bures metric 
d{-, ■) 2^ that, for pure states, is given by 



d($,*)= min (5($,e'"*) = ^2 - 2|(4>, *)|. (8) 

C(e[0,27r] 

Also, we need to compare distributions, and we consider Hellinger metric 23| and a natural gen- 
eralization to several distributions. In Appendix B, explicit expressions of these metrics and their 
properties can be found. 

The physical imposition operator needs, as input, sets of compatible distributions (in order to 
agree with the indeterminacy principle), and in computer simulations we need a way to obtain 
them. This is why it is important to define the notion of generator state: Let $ £ "H be a quantum 
state and {A^ }j=i..m be a set of observables with eigenvectors bases {<p{.}fc=i..Ar respectively. Then, 
we are able to generate the set of distributions 

Pi, = |(^^„$)|2,fc-l..iV,j = l..m, (9) 

from the state $. After this process, we remember the distributions p'f. but we forget the generator 
state $, and we try to reconstruct it using the physical imposition operator, the eigenvectors 
bases, and the information contained in the distributions {p^}. Generator state appears as an 
useful subscript in the physical imposition operator Tab$, but we remark that the state $ is not 
a known state in the reconstruction process. We only have knowledge about the distributions 
generated from it. 



B. Quantum state reconstruction from dynamical system 



In this section, we show that the physical imposition operator has a beautiful description from 
dynamical system theory, and we prove that each Pauli partner is an attractive fixed point of the 
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physical imposition operator. First, let us introduce some basic notions from dynamical systems 
theory 



24 



25|. 



DEFINITION III.2 (Fixed point) Let cf) e U and T : % ^ % he an operator. We say that (f> 
is a fixed point ofT iff (j) is invariant when T is applied. That is, T(j) — (p. 

DEFINITION III.3 (Attractive fixed point) LetT -.U^U, d{-, ■) a metric defined on the 
rays of % and cf) a fixed point of T . We say that (p is an attractive fixed point of T if d{Til), (f) < 
d{ip,(j)) for all -0 contained in a neighborhood of (j). 

It is easy to notice that every fixed point of a single physical imposition operator has the same 
distribution than the generator state. Then, a state is a fixed point of the single physical imposition 
operator iff it is a Pauli partner for the single observable considered. However, if we consider more 
than one observable in the physical imposition operator there are more fixed points than partners. 
That is, there could exist a state 77 satisfying Tas*?? — TB<i>TA<i>ri — rj, but Ta$?7 ^ rj (notice that 
TAB<i>fl — 1] necessarily imply that Ts'^rj = rj, because Tg* is the last single operator applied and it 
impose the eigenvalues distribution about B contained in $). Let us define the set of fixed points 
that are interesting for the quantum state reconstruction problem. 

DEFINITION III. 4 Let A be an observable and $ a generator .state. We define as the set 
of fixed points of the single imposition physical operator . That is, 

F^i = e •H/Ta** = *}. (10) 

Here, we are considering only one representant 5* from the ray — e'"^. Now, we are going to 
define a particular set of fixed point of the multiple physical imposition operator. 

DEFINITION III. 5 Let {A^}j^i,,m be a set of incompatible observables and <^ £% a generator 
state. We define F^i ^771 ^ as 



n-.-nF^™*. (11) 



Notice that since $ is a generator of distributions we always have IF^i.-.a™*! 7^ where the symbol 
I • I is the cardinality of the set. The next two propositions are easy to understand and their purpose 
is to clarify the recent definitions. The proofs are trivial. 

PROPOSITION III.l Let {A-'jj^i .m be an informationally incomplete set of observables and 
$ G H a generator state. Then, the number of Pauli partners M is given by 

AA=|Fai...a."*I- (12) 
In case of Af = I we say that $ is a Pauli unique. 

PROPOSITION III. 2 A set of observables {^-'}j=i..m is informationally complete iff 

|Fai...a"*I = 1, V$eH. (13) 
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PROPOSITION III. 3 Let {^-'}j=i..m be a set of m-observables and $ G "H a generator state. 
Then, Xl/ G r^i...^m$ 'ijf ^ is a Pauli partner of ^. 

Now we are going to present the most important fact on this paper: 

PROPOSITION III. 4 Let {A^}j^i, „i be a set of m-observables, $ G "H a generator state and 
"i! G r^i...^77i^. Let •) be the Bures metric for quantum states and r4i...^m$ the physical 
imposition operator related to the observables {^-'}j=i..m- Then, 5* is an attractive fixed point of 

The proof of this proposition is easy but not short, and it can be found in Appendix A. This 
proposition means that all fixed points in Ty^i j^m^, that is, the complete set of Pauli partners of 
the system, are attractive fixed points of T^i considering Bures Metric. The most important 
result given by the above proposition is that the multiple physical imposition operator allows 
us to find the complete set of Pauli partners of a system for every set of observables and every 
generator chosen. The above proposition together with the fact that probability distributions 
come from a generator state are necessary and sufficient conditions to have a convergent sequence 
= T^i ^m$4'o- In Section IV we present some results obtained from numerical simulations. 
We can extend the definition of the physical imposition operator to mixed states and compos- 
ited systems, considering tensor product of singles ones. Another interesting case is to consider 
entangled physical imposition operators, entangled generators or both of them. A further study of 
our algorithm in these cases will be presented in another work. 

C. Bifurcations 

Sometimes, for some particular values of the parameters of an operator, the stability of their 
fixed points change, fixed points are added or removed. When one of these situations happen we 
say that the operator has a bifurcation 24] . In our work, bifurcations are very important, because 
they tell us about where Pauli partners appear or disappear. The parameters in the physical 
imposition operator are given by the physical distributions. Let us give an example of bifurcations 
in the physical imposition operator: let {A, B) be two observables with MUBs eigenvector bases, 
defined in a N dimensional Hilbert space where iV is a power of a prime number. Let us consider 
the following generators: $1 be an eigenvector of A or B, and $2 a mutually unbiased vector to 
A and B eigenvectors. This means that $1 generates (sharp, flat) or (fiat, sharp) distributions on 
{A,B) whereas $2 generates (fiat, flat) distributions on {A,B). We know that $1 has no partners 
for the considered observables, because it is an eigenvector of A or B, whereas $2 has at least 
N{N — 1) — 1 partners corresponding, together with $2, to the set of eigenvectors of the — 1 
complementary MUBs to the eigenvectors of A and B. Then, every curve in the state space 
connecting the generators $1 with $2 contains at least one bifurcation. 

The basin of attraction of Pauli partners has an interesting property. Let us notice that the 
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sequences 

*n = r;BC...$*0, (14) 

and 

*n = nBC...$*0, (15) 

are identical when T/i$^'o = T^^^'o, that is, when ^Pq have the same set of unitary phases than ^'o 
in the basis of eigenvectors of A, namely {<Pfc}fe=i..Ar. This means that the information about the 
basin of attraction of Pauli partners is only contained in the — 1 relevant phases of the initial 
state ^I'o- Then, if we consider the decomposition of $o in the basis {(fik} 

*o = (^/^,VP^e'"^•••,VP^e^""")> (16) 
the — 1 dimensional real vector 

f = (ai, • • • ,aAr - 1), (17) 

contains the complete information about where the sequence = ^abc -^o converges. This fact 
is an advantage in numerical simulations, because we should consider a — 1 dimensional seed r 
instead of the 2{N — 1) dimensional seed ^'o- 

Until now, we have been studying Pauli partners, that is, considering informationally incomplete 
sets of observables. Let us present a proposition that connects bifurcations with informationally 
complete set of observables in a univocal way. 

PROPOSITION III.5 LetA^..A"' be a set of incompatible observables. Then, Tak has no 
bifurcations for all ^ (E H iff {A^ ..A'^} is an informationally complete set of observables. 

Proof: Suppose that T^i j^m^ has no bifurcations. We know that every eigenvector of A^..A^ has 
no Pauli partner, because there is a unique state that has a sharp distribution. Since T^i y^m^ has 
no bifurcations, Pauli partners cannot appear for any generator state and the quantum state can be 
reconstructed efficiently for all distributions. Then, A^..A™' is an informationally complete set of 
observables. The reciprocal argument can be immediately proved by definition of informationally 
complete set of observables. I 
Notice that this proposition imply a high cost in numerical simulations because we should con- 
sider a huge set of generators. However, some theorems related to necessary conditions for the 
existence of bifurcations could be very useful 1261 ] . Interesting theorems related with information- 



ally complete sets can be found in [9|, |27h2' 
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IV. NUMERICAL RESULTS 



A. Mutually unbiased bases (MUBs) 



We know that it is possible to construct a maximal set of + 1 MUBs in dimension N = p'^, 
where p is a prime number and r is a positive integer. For non power of prime dimensions the 
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maximal set of MUBs is still unknown. For example, in dimension iV = 6 three MUBs have 
been found, but not 7, and in general we know that there exists at least 3 MUBs for every 
dimension. This last case is particularly interesting in physics because, the 3 MUBs are conformed 
by the eigenvector basis of the shift generator = e^^^^ , the boost Bp — e^^-^ and the operator 
ZxBp = p-ixoP ^ipoX _i_ ^i{poX-xoP) (iggjied in a finite dimensional Hilbert space of dimension N 
[soi l. Notice that the last equality is not allowed in finite dimensions, but it is valid when N ^ oo. 
The problem of determining the maximal number of MUBs is contained in the quantum state 
reconstruction problem. For example, given two observables with mutually unbiased eigenvectors 
basis, like position and momentum, we can ask to our algorithm about the maximal set of states 
having fiat distributions for these observables. If TV is prime, the algorithm presented in this work 
is an useful tool to reconstruct the complete set of + 1 MUBs. Starting with position and 
momentum eigenvectors bases we could found + 1 MUBs for = 2, . . . , 37 {N a prime number). 
The algorithm is very efficient. For example, we could reconstruct 24 MUBs in = 23 in 27 
seconds and 38 MUBs in = 37 in 75 minutes in a standard Pentium IV. When A^ is not prime, 
starting with position and momentum basis we could only obtain 3 MUBs. This means that the 
second basis considered (momentum basis) is not adequate for constructing A^ + 1 MUBs for all 



dimensions (if they exist) [3l[. We know that if N = (power of prime) a suitable second basis 
is 4>p = ^]=i where {4ip}p=i..N is the momentum basis in dimension p and is the tensorial 
productory. This second basis is a very particular option, and it cannot be found from searches 
realized over random orthogonal bases. Of course, when we put the correct second basis , the 
algorithm is able to reconstruct A^+1 MUBs when A^ is a power of a prime number, but this process 
is not possible without the additional information about the second basis. Several simulations had 
been realized for A^ = 6 and we could not find more than 3 MUBs, but these numerical results are 
not enough to refuse the existence of 7 MUBs in A^ = 6. The problem of determining maximal sets 
of MUBs in non power of prime dimensions starts in the very huge set of options for the second 
basis. Moreover, the choice of the first vector of the second MUB is determinant to construct a 
maximal set of MUBs. Considering the first MUB as the canonical (computational) basis, every 
vector of the second basis has A^ — 1 complex phases to be fixed, because all the amplitudes are 
equal to 1/\/N. The additional condition needed is the orthogonality of the bases — N/2 

conditions). Then, given that we have A^(A^ — 1) phases, there exist a ~ A^/2-manifold for 

the free parameters, namely {xk}k=i..N^/2-N/2- In numerical simulations, we must take a finite 
partition of each parameter Xk in order to explore all possible second basis. Considering that every 
partition is regular and it has m intervals, then we should analyze /2~n/2 possibilities for the 
second MUB. When A^ = 6, and considering m = 5, which is not a good partition, we should 
analyze 5^^ w 3 x lO^'' possibilities. Supposing that wc have an algorithm that decides if 2 MUBs 
can be extended to A^ + 1 MUBs in only 1 second we need wait 968 thousand years to analyze every 
case. Because of this problem, numerical existence of maximal sets of MUBs should be explored 



in a different way, for example, considering MU constellations 
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B. Basin of attraction 

Finally, we present some numerical results for the basin of attraction of the physical imposition 
operator. In simulations we could see that basin of attractions are very complex sets in general, 
and due to this complexity we do not have any hope of obtaining in a simple way analytical 
expressions neither for the basin of attraction of Pauli partner nor for the undesirables initial 
states. Undesirable initial states are initial states 4*0 within the basin of attraction of an attractive 
fixed point that is not a Pauli partner. In our several simulations realized in dimension iV = 3 we 
could not detect basins of attraction with fractal behavior, and we have no knowledge about the 
existence of fractal basin of attraction in higher dimensions. 

As we have seen in Section IIII C[ we can show in a plane the complete basins of attraction for 
every set of observables and every generator when iV = 3. Let us consider, then, two observables 
with MUBs eigenvectors basis defined on H = C^, specifically position and momentum. In Fig. 1 
we show the basin of attraction of Pauli partners in case of the MUBs problem, that is, considering 
a generator of flat distributions. We found 6 partners, corresponding to the two complementary 
mutually unbiased bases to the eigenvectors bases of the observables. This result tells us that 
there is only one way to complete a maximally set of MUBs from position and momentum when 
= 3. Another interesting result is that there is no isolated MU partner, though they exist in 
other lower dimensions [s^. Fig. 2 and Fig. 3 correspond to two random generators and we find 
4 and 5 Pauli partners respectively. In Fig. 4 and Fig. 5 we consider two random generators 
without partners. Each color in the figures is identifying a different basin of attraction. Violet 
regions shows basin of attraction of the undesirable fixed points. Sometimes, an undesirable point 
may be a saddle point instead of a fixed point, having an attractive manifold of null measure, 
named stable manifold. Applying our algorithm to the maximally MUBs problem (N = 3) we 



could see that every basin of attraction of a partner is a triangle of identical area [38|. In Fig. 
1, undesirable states are contained in a null measure stable manifold living in the border of the 
triangles. This null measure attractive manifold can be numerically detected due to the nature 
of the convergence process in numerical simulations, that is, a null measure manifold becomes a 
little set with the same dimension than the Hilbert space, and the size of these sets is a growing 
function of the numerical convergence bound condition. Unfortunately, undesirables states in the 
MUBs problem, in case of A > 3, are fixed points and consequently they have a basin of attraction 
that has not null measure. Moreover, the number of undesirable fixed points increases faster than 
the number of Pauli partners with the dimension N. For example, in case of A^ = 31 most of time 
the algorithm is rejecting undesirable states. This is the reason why we do not have a polynomial 
algorithm for quantum state reconstruction. If undesirable states could be avoided the algorithm 
would be polynomial (linear in case of MUBs problem), and the quantum state reconstruction 
problem would be reduced to a NP-EASY problem. 
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V. CONCLUSIONS 



We have presented an algorithm for state reconstruction in finite dimensional Hilbert spaces. We 
showed that every Pauli partner is an attractive fixed point of the physical imposition operator, 
independently of the set of observables considered. This fact allows us to obtain the complete set of 
Pauli partners for a given system. The algorithm was applied to the following particular cases: 1 ) 
the MUBs problem, where we could reconstruct the maximal set of MUBs from N = 2 to N = 37 
{N prime), 2) the study of basin of attraction of Pauli partners in case of A'' = 3 considering 
position and momentum observables. An interesting property found for the physical imposition 
operator is when it has no bifurcations, because this is a necessary and sufficient condition to have 
an informationally complete set of observables. Unfortunately, our algorithm can not be applied, 
at least in the way considered, to contribute to numerical evidences of the non existence of a 
maximal set of MUBs in non power of prime dimensions, and this is due to the high dimensional 
manifold that characterizes the different options for the second MUB. The physical imposition 
operator for mixed states and composited systems will be considered in another work. Extensions 
of our algorithm will be studied in order to try to find numerical solutions to the problems of MU 
Constellations, SIC-POVMs, Equiangular Tight Frames and Quantum t-Designs. 

We thank to Dr. T. Santhanam for sending us his works about quantum mechanics in finite 
dimensions. This work contains the last results of the PhD Thesis of DG, supported by a CONICET 
scholarship. This work was supported by CONICET and CONICyT PFB-0824. 
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Appendix A 



Our intention here is to give a proof of Proposition IIII.41 First, we need to proof some previous 
geometrical properties for the single physical imposition operator: 

PROPOSITION A.l Let A he an observable, $ G H a generator state, Ta^ the physical impo- 
sition operator and d{-,-) the Bures metric. Then, 

1. ^(Ta**,*) < V$,* e n. 

2. d{TA^-^,Vk) = d{^,^k). Vfc = l...A^. 



3. d(TA<i>*,$) < 2minfcd($,(^fc) = 2V2^1 - max/, VI' G 

I diTA^-^, $) < 2d(*, $), V $, e 

5. Ta^^ = e ^ < V ^' e iVA(C), 

where Na{S,) is a neighborhood of ^. 
proof: 



1. 



N 



fc=l 



fc=l 

AT-l 



> 



Then, 



< V2Vl-|(v|/,(f>)| 

2. Remembering the definition of the physical imposition operator we can find that 



\{TA^'^,ipk)\ - 



= l(<f,¥'fe) 



Then, 



d{TA^^,^k) - V2v/l-|(TA**,<Pfc)| 
= d{^,ipk) Vfc = l...A^. 



(Al) 

(A2) 

(A3) 
(A4) 

(A5) 
(A6) 
(A7) 

(A8) 
(A9) 

(AlO) 
(All) 
(A12) 
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3. Using triangle inequality and Ea. (IA12p 

d(T^**,$)<d(T^$*,^fe) + d((^fe,$) = 2d(^fe,$), yk^l...N. (A13) 
The more restrictive condition is given by 

d(rA<i.*,$) < 2mind((^fc,$). (A14) 

k 

Equation 

niind(<^fc,$) = V2./I - max^/p^, 

k \ k 

is proven immediately from Bures metric definition. 

4. From triangle inequality and property 1 ) from this proposition we have 

ci(rA$«',$) < ^(Ta**,*) + (A15) 
< 2d(*,$). (A16) 

5. Taking into account Ea. (|A7p and doing the following parameter change ^ -> ^ + (5^, $ — ^ 
we have 

rf(T^$(C + 50, ^ + 50 < + St 0- (A17) 

Notice that in the last equation we have written instead of Ta^- This is in order to 
consider the general situation, when $ is a partner of ^ and not necessarily Taking 

\{Ta^{^ + 50,^ + S0\ = \{TA<!>{i + so,0 + {TA^i>{i + so,so\ 
« \{Ta<!>{^ + so,01 

or, equivalently 

d{TA<s.{^ + S0,^ + SO - diTA^ii + 50, 0- (A18) 
Taking into account Eas. (|A17p and (|A18p we have 

d(TA$(e + 5c),e) « d(r^*(e + 5e),e + 5c) (ai9) 

< d{^ + 600- (A20) 



Item 4 ) in previous proposition defines a relationship between the distance between two elements 
before and after physical imposition operator is applied. The factor 2 on this equation is related 
with the existence of Pauli partners. A factor less than one would mean a contradiction to existence 
of Pauli partner, because, in this case, the physical imposition operator would be a contraction 
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and, by Banach's fixed point theorem, it would have a unique fixed point, rejecting the idea of 
Pauli partners. 

Now, we are able to give a proof of ProDOsition lIIL4l Considering the last proposition for several 
observables A,B,C, - ■ ■ we have 

d{^,0 > d{TA^^,0 (A21) 

> d(TB$T^*^',0 (A22) 

> d(Tc*TB*r^$*,0 (A23) 
: (A24) 

> d{---Tc<!>TB^TA^^,0 (A25) 

> d(TABC...**,0, (A26) 

where we consider that * G Nabc---{0 = Na{£.) n Nb{0 ^ ^ciO ■■■ ■ Then, £_ is an attractive 
fixed point of Tabc - '^- Notice that the neighborhood Nabc---{0 cannot be an empty set or a set 
of null measure in state space, because each set contain an open set around ^. Given that an 
intersection of a finite number of open sets is an open set, then Nabc---{0 contain, at least, an 
open set. Notice that the basin of attraction of the multiple physical imposition operator contain 
the set Nabc-{0- ■ 



Appendix B 

In this appendix we will define two useful metrics in quantum mechanics. Suppose that we have 
two different states $, e "H and we want to know if they have the same distribution for the 
eigenvalues of some observable A. Several metrics can be defined for doing that, and a nice option 
is Hellinger metric !?(•,•). 

DEFINITION B.l Let A be an observable defined on a N -dimensional Hilbert space %, 
{'fk}k=i..N its eigenvectors basis and $,'1' G Then, Hellinger metric is given by 

/ N \ 1/2 

DA{<i>,^)= [YSVf^-V^fj , (Bl) 
where pk = |((/?/c,$)P and ak = |((^/c,*)P- 

Hellinger metric is able to compare distributions for only one observable. However, Pauli partners 
has the same distributions for several observables and we need to define a metric that considers all 
of them. In order to do this, let us define an intuitive generalization of Hellinger metric, that we 
named distributional metric. 

DEFINITION B.2 Let A^..A"^ be a set of m-observables defined on a N -dimensional Hilbert 
space H, {'/'feli^i"'^ their eigenvectors basis and $,4' G Then, the distributional metric is 
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given by the expression 

/ m \ 
Va.„a^{^,^)= I^P^{Da,{<^,^))A . (B2) 

It is easy to prove that the distributional metric satisfies all metric conditions. Moreover, the 
distributional metric is proportional to the usual metric in . Let us establish an upper bound 
for distributional metric considering Bures metric. 

PROPOSITION B.l Let be a set of m-observables and two different states. 

Then, Bures metric is an upper bound for every Bellinger metric and for the distributional metric. 
That is, 

> Dao ($, V * e ^, Vi = l..m (B3) 

and 

> ©Ai.. A- V$,* e ^,Vyl^.A™. (B4) 

Proof: 

Let A be an observable with eigenvector basis {vfe}fe=i..jv, The decomposition of the 

states $ and ^ in the eigenvector basis ^pk is 

N 

* = l]v^e^"''¥'fe (B5) 

fc=i 

N 

* = l^v^e^^'-y'fc- (B6) 

fe=i 

Using the triangular inequality we can see that 



= |^^/p^^/^e^('3^-«'=)| (B7) 
fe=i 

N 

< I^|\/A^x/^e^('3'=-«'=)| (B8) 
fe=i 

N 



fe=i 

Remembering the expresion of the Bures metric 



= \/2Vl- K$,*)|, (BIO) 

and Hellinger metric 



N 



1 - I] V^V^, (Bll) 
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we obtain 

> £1^(4', V$, * e -H. (B12) 

Then, Eq.(|B3| is proven. 

Now we are going to proof Ea. (jB4[) . The summatory of the m-inequahties Eq. (jB3l) considering 
every m-observable is given by 

m 

md^<i>,^)>J2iDA.{'^,^)f, (B13) 



or, equivalently 



\ 



-^(^?^4'i>,*))^ (B14) 

T7 ^- ^ 



m 



Then, 



> 2?^i..An.($,^'), V$,* e -H,VA\.A'". (B15) 

I 

This proposition is another manifestation of the existence of Pauh partners, because 

= ^ i:>Ai..A-($,*) = 0, (B16) 

but 

I^Ai^A". ^ = 0. (B17) 

However, Ea. (jB17|) becomes an imphcation when A^..A™ is an informationally complete set of 
observables. For example, it is valid when we consider a complete set of + 1 mutually unbiased 
observables in an N dimensional Hilbert space (iV-power of prime). 
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FIGURE CAPTIONS (Figures attached) 

FIG 1: Basin of attraction of Pauli partners in case of position and momentum observables with 
a flat generator {N = 3). The six partners found form two orthogonal MU bases. 

FIG 2: Basin of attraction of Pauh partners in case of position and momentum observables with 
a random generator state (4 partners found). Violet regions correspond to basin of attraction of 
undesirable fixed points. 

FIG 3: Basin of attraction of Pauli partners in case of position and momentum observables with 
a random generator state (5 partners found). Violet regions correspond to basin of attraction of 
undesirable fixed points. 

FIG 4: Basin of attraction of Pauli partners in case of position and momentum observables with 
a random generator state. In this case there is a Pauli unique state. Violet regions correspond to 
basin of attraction of undesirable fixed points. 

FIG 5: Basin of attraction of Pauli partners in case of position and momentum observables with 
a random generator state. In this case there is a Pauli unique state. Violet regions correspond to 
basin of attraction of undesirable fixed points. 
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